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Abstract 

This  paper  reports  a  nonlinear  modeling  study  of  a  solid  oxide  fuel  cell  (SOFC)  stack  using  a  least  squares  support  vector  machine  (LS-SVM). 
SOFC  is  a  nonlinear,  multi-input  and  multi-output  system  that  is  hard  to  model  by  traditional  methodologies.  So  far,  most  of  the  existing  models 
are  based  on  conversion  laws,  which  are  very  useful  for  cell  design.  However,  they  are  too  complicated  to  be  applied  to  control  system  design. 
To  facilitate  a  valid  control  strategy  design,  this  paper  tries  to  avoid  the  internal  complexities  and  presents  a  black-box  model  of  the  SOFC  based 
on  LS-SVM.  The  simulation  tests  reveal  that  it  is  feasible  to  establish  the  model  using  LS-SVM.  At  the  same  time,  the  experimental  comparisons 
between  the  LS-SVM  model  and  radial  basis  function  neural  network  (RBFNN)  model  demonstrate  that  the  LS-SVM  is  superior  to  the  conventional 
RBFNN  in  predicting  stack  voltage  with  different  fuel  utilizations.  Furthermore,  based  on  this  black-box  LS-SVM  model,  valid  control  strategy 
studies  such  as  predictive  control,  robust  control  can  be  developed. 

©  2006  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

The  solid  oxide  fuel  cell  (SOFC)  is  an  energy  conversion 
device  that  produces  electricity  by  electrochemically  combin¬ 
ing  fuel  (e.g.  hydrogen)  with  oxidant  (e.g.  oxygen)  gases  across 
an  ionic  conducting  oxide  [1].  SOFC  provides  many  advantages 
over  traditional  energy  conversion  systems  including  high  effi¬ 
ciency,  modularity,  fuel  adaptability  and  very  low  levels  of  NOx 
emission.  An  important  tool  in  fuel  cell  development  is  mathe¬ 
matical  modeling,  which  is  particularly  appropriate  for  SOFCs, 
where  localized  experimental  measurements  are  difficult  due  to 
the  high  operating  temperature  [2].  The  results  obtained  from  a 
reliable  and  effective  model  can  be  very  useful  to  guide  future 
research  for  fuel  cell  improvements  and  optimization. 

It  is  well  known  that  the  SOFC  system  is  sealed,  and  works 
in  a  complicated  high-temperature  (600-1000  °C)  environment. 
As  a  nonlinear  multi-input  and  multi-output  system,  SOFC  is 
hard  to  model  using  traditional  methodologies.  In  the  last  sev¬ 
eral  decades,  fruitful  results  from  SOFC  stack  modeling  have 
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been  obtained  [3-5].  However,  most  of  the  existing  models  that 
have  been  developed  emphasized  the  detailed  description  of  cell 
internal  processes,  such  as  mass  balances,  energy  balances  and 
electrochemical  kinetics.  These  models  are  very  useful  for  cell 
design,  but  they  are  too  complicated  to  be  used  for  a  control 
system  design. 

To  meet  the  demands  of  developing  valid  control  strategies, 
some  researchers  have  attempted  to  establish  novel  SOFC  mod¬ 
els.  A  black-box  identification  technique  such  as  the  artificial 
neural  network  (ANN)  has  been  used  to  derive  a  SOFC  model 
from  the  experimental  data  quickly  [6].  Although  this  ANN 
model  shows  a  high  accuracy  and  is  much  faster  and  easier  to 
use,  its  practical  design  suffers  from  drawbacks  such  as  the  exis¬ 
tence  of  local  minima  and  over- fitting,  choice  of  the  number  of 
hidden  units,  etc.  So  a  new  modeling  approach  is  needed  to  pro¬ 
vide  a  better  solution.  In  this  work,  a  least  squares  support  vector 
machine  (LS-SVM)  is  presented  to  establish  a  black-box  model 
for  the  SOFC. 

LS-SVM  proposed  by  Suykens  and  Vandewalle  [7]  is  a  mod¬ 
ification  of  the  standard  SVM.  Unlike  ANN,  LS-SVM  possesses 
prominent  advantages:  over- fitting  is  unlikely  to  occur  by  adopt¬ 
ing  the  structural  risk  minimization  (SRM)  principle,  and  the 
global  optimal  solution  can  be  uniquely  obtained  by  solving  a 
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Nomenclature 

b 

bias  term 

e 

predictive  error 

/(•) 

nonlinear  function 

/ 

stack  current 

K(x,  x{) 

kernel  function 

m 

lag  of  the  fuel  utilization 

n 

lag  of  the  output  voltage 

N0 

number  of  cells  in  the  stack 

P 

stack  dc  output  power 

<7h2 

input  hydrogen  flow  rate 

6102 

input  oxygen  flow  rate 

T 

stack  operating  temperature 

u 

fuel  utilization 

V 

stack  output  voltage 

V 

predictive  output  voltage 

w 

weight  vector 

U 

ith  original  training  datum 

•Unax 

maximum  of  the  training  data 

Unin 

minimum  of  the  training  data 

x' 

normalized  result 

Greek  letters 

Oli 

Lagrange  multipliers 

y 

regularization  parameter 

a 

kernel  width 

set  of  linear  equations.  A  number  of  structures  and  algorithms  for 
modeling  using  LS-SVM  have  been  proposed  [8-10].  However, 
the  concrete  study  of  modeling  SOFC  with  LS-SVM  cannot  be 
found  in  prior  papers. 

This  paper  is  organized  as  follows.  In  Section  2,  a  brief  analy¬ 
sis  of  the  characteristics  of  a  SOFC  stack  is  presented.  In  Section 
3,  LS-SVM  for  nonlinear  system  modeling  is  explained.  In  Sec¬ 
tion  4,  Identification  structure  of  a  SOFC  stack  and  the  detailed 
processes  of  training  and  testing  the  LS-SVM  model  are  given. 
In  Section  5,  conclusions  and  suggestions  for  future  work  are 
summarized. 

2.  Description  and  analysis  of  SOFC  stack 


Fuel  Cell  Stack 


Anode 

Flow 


Cathode 

Flow 


Fig.  1.  Single  cell  and  stack  of  SOFC. 

an  external  circuit  [11].  The  typical  hydrogen  fed  cross-flow 
configuration  of  a  SOFC  is  shown  in  Fig.  1  [12]. 

A  single  cell  produces  an  open-circuit  voltage  of  approx¬ 
imately  1  V.  Cells  have  to  be  connected  together  in  a  series 
arrangement  to  form  a  cell  stack  that  delivers  the  higher  voltages 
suited  to  static  converters  [13]. 

As  we  know,  cell  voltage  calculation  is  the  core  of  any  fuel 
cell  modeling.  For  a  given  SOFC  stack,  the  output  voltage  is 
influenced  by  many  operating  parameters  such  as  temperature, 
pressure,  fuel  utilization,  flow  rate,  etc.  However,  due  to  the 
high  number  of  operating  variables,  a  complete  experimental 
database  of  SOFC  under  the  different  operating  conditions  is 
difficult  to  obtain  and  no  data  are  available  in  the  open  litera¬ 
ture  yet  [14].  Up  to  now,  almost  no  model  has  ever  been  able  to 
accommodate  all  these  operating  variables.  Our  LS-SVM  model 
is  no  exception.  Fuel  utilization  is  one  of  the  most  important 
operating  parameters  for  a  fuel  cell  and  has  significant  effects 
on  the  cell  voltage.  In  order  to  analyze  the  effects  of  different 
fuel  utilizations  on  output  voltage,  we  chose  current,  which  is 
determined  by  the  external  load,  and  fuel  utilization  as  variables, 
while  holding  other  operating  parameters  constant.  Based  on  the 
LS-SVM  approach,  we  present  a  voltage-fuel  utilization  model 
under  different  currents  in  this  paper.  Furthermore,  like  the  stan¬ 
dard  S VM,  LS-SVM  also  has  better  generalization  performance 
and  this  ability  is  independent  of  the  dimensionality  of  the  input 
data.  So  our  LS-SVM  model,  obtained  with  the  two  variables, 
can  predict  stack  voltage  as  precisely  as  a  model  considering 
more  variables.  Besides,  by  adding  more  variables  to  our  LS- 
SVM  model  and  training  it  again,  the  new  multi- dimensional 
model  can  be  obtained  easily. 


A  brief  description  and  analysis  of  the  SOFC  stack  is  given 
in  this  section.  The  reason  for  establishing  a  cell  voltage-fuel 
utilization  model  under  different  cell  currents  is  also  explained. 

A  SOFC  consists  of  an  interconnected  structure  and  a 
three-layer  region  composed  of  two  ceramic  electrodes,  anode 
and  cathode,  separated  by  a  dense  ceramic  electrolyte  (often 
referred  to  as  the  PEN,  Positive-electrode/Electrolyte/Negative- 
electrode).  In  this  cell,  the  oxygen  ions  formed  at  the  cathode 
migrate  through  the  ion-conducting  solid  ceramic  electrolyte  to 
the  anode/electrolyte  interface  where  they  react  with  the  hydro¬ 
gen  and  carbon  monoxide  contained  in  (and/or  produced  by) 
the  fuel,  producing  water  and  carbon  dioxide  while  liberating 
electrons  that  flow  back  to  the  cathode/electrolyte  interface  via 


3.  LS-SVM  for  nonlinear  system  modeling 

In  the  following,  we  briefly  introduce  LS-SVM  algorithm  for 
nonlinear  system  modeling,  based  on  [7,15]. 

Assume  a  set  of  training  data  is  given: 

(xi,yi),...,(xN,yN)£R"xR  (1) 

The  nonlinear  function  ^(-)  is  employed  to  map  the  orig¬ 
inal  input  space  Rn  to  high  dimensional  feature  space 
=  (cp(x i),  (p{x 2), . . .,  cp(x n)).  Then  the  linear  decision  func¬ 
tion  y(xt)  =  wT(p(xi)  +  b  is  constructed  in  this  high  dimensional 
feature  space.  Thus  nonlinear  function  estimation  in  original 
space  becomes  linear  function  estimation  in  feature  space. 
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The  quadratic  loss  function  is  selected  in  LS-SVM.  Then  the 
optimization  problem  of  LS-SVM  is  formulated  as: 

1  1  N 

min  J(w,  e)  =  -wTw  +  y~y^ej,  y>  0  (2) 

w,b,e  2  2Z — ' 

i= 1 


The  selection  of  kernel  function  needs  some  knowledge  in 
advance,  there  is  no  common  conclusion  currently.  In  this  paper, 
the  RBF  function  is  used  as  the  kernel  function  of  LS-SVM 
because  RBF  kernels  tend  to  give  good  performance  under  gen¬ 
eral  smoothness  assumptions. 


subject  to  the  equality  constraints 

yi  =  wT(p(xi )  +  b  +  ei,  i  —  1, N  (3) 

We  construct  the  Lagrangian  as 

N 

L(w,  b,  e,  a)  =  J(w,  e)  -  y^afiw1  (p{xj)  +  b  +  et  -  yt}  (4) 


i= 1 


where  oti  (i  =  1, . . .,  AO  are  the  Lagrange  multipliers.  The  condi¬ 
tions  for  optimality  are  given  by 


dL 


=  0- 

>  w  = 

dw 

dL 

=  0  - 

>  oti  = 

det 

dL 

=  0- 

->  yi  = 

N 


i= 1 


dL 

~db 


N 


=  o  ->  =  o, 


i=l 


da 


(5) 


With  solution 


'0  1  1 

b 

1  K(x\,  X\)  +  \/y  ■■■  K{x \,XN) 

...  • 

a\ 

•  •  *  • 

.  1  K(xn,x  1)  •••  K(xN,xN)  +  l/y_ 

.OiN  _ 

'  0  ' 

_yN  _ 


The  resulting  LS-SVM  model  for  nonlinear  system  becomes 

N 

y(x)  =  ffdiKix,  Xi)  +  b  (7) 

i=  1 

where  au  b  are  the  solution  to  the  linear  system.  Using  the  nor¬ 
mal  linear  equations  program  method,  we  can  get  the  parameters 
a.i  of  (6).  By  the  Karush-Kuhn-Tucker  (KKT)  conditions,  the 
parameter  b  can  be  calculated,  so  the  LS-SVM  model  for  non¬ 
linear  system  can  be  obtained. 

The  kernel  function  K(x ,  x{)  is  any  symmetric  function  that 
satisfies  Mercer’s  condition.  The  typical  examples  of  kernel 
function  include  linear,  polynomial,  radial  basis  function  (RBF) 
kernel. 

T 

Linear  :  K(x i,  X2)  =  x[X2  (8) 

Polynomial  :  K(x i,  X2 )  =  ( xjx2  +  l)P ,  p  e  N  (9) 


4.  Modeling  SOFC  based  on  LS-SVM 

A  LS-SVM  can  be  regarded  as  a  black-box  which  can  pro¬ 
duce  certain  output  data  as  a  response  to  the  specific  input  data. 
In  this  modeling  procedure,  the  relationship  between  input  and 
output  of  SOFC  can  be  emphasized  while  the  sophisticated  inner 
structure  is  ignored.  In  order  to  establish  the  expected  nonlinear 
model  of  SOFC,  we  choose  fuel  utilization  and  cell  current  as 
the  model  inputs,  and  cell  voltage  as  the  output.  In  the  follow¬ 
ing,  identification  structure  of  SOFC  stack  based  on  LS-SVM  is 
given  firstly,  and  then  the  processes  of  training  and  testing  the 
LS-SVM  model  are  presented. 

4.1.  Identification  structure  of  SOFC  stack  based  on 
LS-SVM 

In  general,  a  wide  class  of  nonlinear  systems  can  be 
described  by  nonlinear  autoregressive  model  with  exogenous 
inputs  (NARX).  So  in  this  paper  the  SOFC  nonlinear  system 
with  two  inputs  and  one  output  can  be  described  as  follows: 

V(k+1  )  =  f[V(k),  V(k-  1),..., 

V(k  —  n),  u(k),  u(k  —  1),  ... ,  u(k  —  m),  I(k)]  (1 1) 

Supposing  there  is  a  series  of  inputs  u(k  —  m),  u(k  —  m 
+  1), . . .,  u(k ),  I(k)  and  outputs  V(k  —  n ),  V(k  —  n  +  1), . . .,  V(k), 
then  the  corresponding  output  V(k+  1)  can  be  obtained  from 
(11).  And  providing  that 

X(k)  =  (V(k),  V(k  —  1),  ... ,  V(k  —  n),  u(k), 

u(k  —  1),  ... ,  u(k  —  m),  I(k)),  k  =  1,  2,  . . . ,  N  (12) 

then 

V(k  +  1)  =  f(X(k))  (13) 

We  firstly  construct  the  training  sample  set  (X(k),  V(k+  1)),  and 
then  the  nonlinear  sample  data  can  be  mapped  as  the  linear  out¬ 
puts  in  high  dimensional  feature  space  by  using  LS-SVM. 
Namely 

N 

V(k  +  1)  =  Y^oiiK(X(k),  X(i))  +  b  (14) 

i= 1 

The  identification  structure  of  SOFC  stack  based  on  LS-SVM 
is  shown  in  Fig.  2,  where  TDL  is  the  tapped  delay  line,  and  the 
predictive  error  e(k  +  1)  =  V(k  1)  —  V(k  1). 

4.2.  Training  process  of  LS-SVM 


RBF  : 


K(x i,  *2)  =  exp 


/-||xi  -x2||2\ 

V  2<r2  ) 


(10) 


In  general,  steps  used  in  training  LS-SVM  include:  train¬ 
ing  data  choosing  and  preprocessing,  selection  of  the  optimal 
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Fig.  2.  Identification  structure  of  SOFC  stack  based  on  LS-SVM. 


LS-SVM  parameters,  testing  data  choosing  and  preprocess¬ 
ing. 

4.2.1.  Training  data  choosing  and  preprocessing 

In  our  study,  a  mathematical  model  in  [16]  is  used  to  generate 
the  data  required  for  the  training  of  the  LS-SVM  model.  The 
mathematical  model  has  been  developed  to  research  the  steady- 
state  feasible  operating  regime  of  the  SOFC.  Here  three  groups 
of  fuel  utilization  and  cell  voltage  data  at  100,  200,  and  300  A 
are  chosen  as  training  data,  and  each  group  has  101  pairs  of  data. 
Main  operational  parameters  of  SOFC  are  varied,  such  as  fuel 
utilization  (0.4-0. 9),  stack  current  (100-300  A)  and  voltage  in 
ranges  that  correspond  to  the  fuel  utilization  and  stack  current 
as  shown  in  Fig.  3.  Some  parameters  of  the  SOFC  stack  used  in 
the  LS-SVM  modeling  are  shown  in  Table  1. 

In  most  cases,  all  given  training  data  are  normalized  to  [0,  1] 
or  [  —  1 ,  1]  in  order  to  increase  the  training  speed,  facilitate  mod¬ 
eling  and  predicting.  In  this  paper,  we  normalize  each  group  of 
training  data,  including  fuel  utilization,  stack  current  and  volt- 


Fig.  3.  Training  data:  stack  voltage  (V)  vs.  fuel  utilization  (u)  at  100,  200  and 
300  A. 


Table  1 

Parameters  of  the  SOFC  stack  used  in  the  LS-SVM  modeling 


Item 

Value 

No 

384 

T 

1273  K 

P 

lOOkW 

I 

100-300  A 

V 

Variable 

u 

0.4-0.9 

<7h2 

1.2e— 3  kmol/s 

<lo2 

2.4e— 3  kmol/s 

age,  to  [0,  1  ]  by 

*'  =  *  Vmin  (15) 

Tmax  -Lriin 

4.2.2.  Selection  of  the  optimal  LS-SVM  parameter 

The  precision  and  convergence  of  LS-SVM  are  affected  by 
regularization  parameter  y  and  kernel  width  o.  So  in  order  to 
obtain  high  level  SOFC  model,  y  and  o  in  the  LS-SVM  have  to 
be  tuned. 

( 1 )  y,  which  determines  the  trade-off  between  minimizing  train¬ 
ing  errors  and  minimizing  model  complexity,  is  important  to 
increase  the  generalization  performance  of  LS-SVM  model. 

(2)  o  influences  directly  the  number  of  initial  eigenval¬ 
ues/eigenvectors.  Small  values  of  o  yield  a  large  number 
of  regressors,  and  eventually  it  can  lead  to  over-fitting.  On 
the  contrary,  a  large  value  of  a  can  lead  to  a  reduced  num¬ 
ber  of  regressors,  making  the  model  more  parsimonious,  but 
eventually  not  so  accurate  [17]. 

Several  researchers  have  presented  some  methods  for  deter¬ 
mining  these  two  parameters,  such  as  bootstrapping,  Bayesian 
methods  and  so  on.  However,  most  of  the  available  methods  can 
be  very  expensive  in  terms  of  computation  time  and/or  training 
data.  For  the  industrial  application  of  LS-SVM,  there  is  a  need 
for  a  fast  and  robust  method  to  estimate  these  two  parameters. 
Fortunately,  we  can  rapidly  tune  these  two  parameters  with  a  10- 
fold  cross-validation  procedure  and  a  grid  search  mechanism  by 
LS-SVM  toolbox  [18].  In  the  final  optimal  LS-SVM  parameters 
are:  y  =  419.0603,  o-  0.4080358. 

4.2.3.  Testing  data  choosing  and  preprocessing 

Testing  data  should  be  different  from  the  data  used  for  train¬ 
ing.  If  testing  data  are  identical  to  training  data,  then  the  LS-SVM 
is  just  interpolating  points  on  a  line — which  is  not  what  we 
expect  the  LS-SVM  to  do  [19].  In  our  study,  the  testing  data 
chosen  for  this  work  are  also  provided  by  the  above-mentioned 
mathematical  model  in  [16].  A  group  of  fuel  utilization  and  stack 
voltage  data  at  280  A  are  chosen  as  testing  data.  Preprocessing 
of  testing  data  is  done  in  the  same  way  as  training  data. 

4.3.  Predicting  with  the  LS-SVM  model 

The  criterion  of  training  LS-SVM  is  to  minimize  sum  squared 
error  (SSE).  After  training,  a  LS-SVM  model  is  obtained,  which 
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Fig.  4.  Voltage-fuel  utilization  characteristics:  predicted  by  LS-S  VM  model  and 
experimental  at  280  A. 


Fuel  utilization 

Fig.  5.  Voltage-fuel  utilization  characteristics:  predicted  by  RBFNN  model  and 
experimental  at  280  A. 

can  be  used  to  predict  new  input  data.  Now  the  trained  LS-SVM 
model  is  used  to  predict  stack  voltage  at  280  A  with  different 
fuel  utilizations.  The  comparison  of  predicted  and  experimental 
voltage-fuel  utilization  curve  at  280  A  is  then  made  to  evaluate 
the  LS-SVM  model’s  prediction  precision  (as  shown  in  Fig.  4). 
At  the  same  time,  RBFNN  model  is  also  used  to  predict  the  stack 
voltage  at  280  A,  and  the  predicted  result  is  shown  in  Fig.  5. 
From  Figs.  4  and  5,  we  can  see  the  LS-SVM  is  superior  to  the 
conventional  RBFNN  in  predicting  stack  voltage  with  different 
fuel  utilizations.  These  indicate  LS-SVM  is  a  powerful  tool  for 
modeling  SOFC  and  our  LS-SVM  model  presented  in  this  paper 
is  accurate  and  valid. 

5.  Conclusions 

To  facilitate  valid  control  strategy  design,  a  nonlinear  model¬ 
ing  study  of  SOFC  using  LS-SVM  is  reported  in  this  paper.  It  is 
shown  that  the  LS-SVM  model  is  more  attractive  in  that  it  avoids 
using  complicated  differential  equations  to  describe  the  stack, 


and  the  input-output  characteristics  can  be  achieved  quickly  by 
LS-SVM  estimation.  Besides,  the  performance  of  our  proposed 
LS-SVM  modeling  approach  has  been  tested  and  compared  with 
the  RBFNN  approach,  and  simulation  results  show  that  the  LS- 
SVM  approach  yields  higher  prediction  accuracy  compared  to 
the  RBFNN  approach.  These  indicate  that  it  is  feasible  to  estab¬ 
lish  the  model  of  the  SOFC  by  using  LS-SVM,  and  the  LS-SVM 
model  presented  in  this  paper  is  accurate  and  valid. 

In  the  future,  based  on  this  black-box  LS-SVM  model,  some 
control  scheme  studies  such  as  predictive  control  and  robust  con¬ 
trol  can  be  developed.  In  addition,  because  better  generalization 
performance  of  LS-S  VM  is  independent  of  the  dimensionality  of 
the  input  data,  a  multi-dimensional  LS-SVM  model  considering 
other  operating  parameters  can  be  obtained  easily. 

Acknowledgment 

This  work  is  supported  by  National  863  Scientific  Project 
Development  Funds  (no.  2003AA5 17020),  PR  China. 

References 

[1]  N.Q.  Minh,  T.  Takahashi,  Science  and  Technology  of  Ceramic  Fuel  Cells, 
Elsevier,  Amsterdam,  the  Netherlands,  1995. 

[2]  P.  Iora,  P.  Aguiar,  C.S.  Adjiman,  N.P  Brandon,  Comparison  of  two  IT 
DIR-SOFC  models:  impact  of  variable  thermodynamic,  physical,  and  flow 
properties,  steady-state  and  dynamic  analysis,  Chem.  Eng.  Sci.  60  (2005) 
2963-2975. 

[3]  E.  Achenbach,  Three-dimensional  and  time  dependent  simulation  of  a  pla¬ 
nar  solid  oxide  fuel  cell  stack,  J.  Power  Sources  49  (1994)  333-348. 

[4]  P.  Aguiar,  C.S .  Adjiman,  N.P.  Brandon,  Anode-supported  intermediate  tem¬ 
perature  direct  internal  reforming  solid  oxide  fuel  cell.  I:  Model-based 
steady-state  performance,  J.  Power  Sources  138  (2004)  120-136. 

[5]  K.P  Recknagle,  R.E.  Williford,  L.A.  Chick,  D.R.  Rector,  M.A.  Khaleel, 
Three-dimensional  thermo-fluid  electrochemical  modeling  of  planar  SOFC 
stacks,  J.  Power  Sources  113  (2003)  109-114. 

[6]  J.  Arriagada,  P.  Olausson,  A.  Selimovic,  Artificial  neural  network  simulator 
for  SOFC  performance  prediction,  J.  Power  Sources  112  (2002)  54-60. 

[7]  J.A.K.  Suykens,  J.  Vandewalle,  Least  squares  support  vector  machine  clas¬ 
sifiers,  Neural  Process.  Lett.  9  (3)  (1999)  293-300. 

[8]  T.  Van  Gestel,  J.A.K.  Suykens,  D.E.  Baestaens,  A.  Lambrechts,  G.  Lanck- 
riet,  B.  Vandaele,  B.D.  Moor,  J.  Vandewalle,  Financial  time  series  pre¬ 
diction  using  least  squares  support  vector  machines  within  the  evidence 
framework,  IEEE  Trans.  Neural  Networks  12  (4)  (2001)  809-821. 

[9]  C.  Lu,  T.  Van  Gestel,  J.A.K.  Suykens,  S.  Van  Huffel,  I.  Vergote,  D.  Timmer¬ 
man,  Preoperative  prediction  of  malignancy  of  ovarian  tumors  using  least 
squares  support  vector  machines,  Artif.  Intell.  Med.  28  (2003)  281-306. 

[10]  B.J.  Kim,  I.K.  Kim,  An  application  of  hybrid  least  squares  support  vector 
machine  to  environmental  process  modeling,  Lect.  Notes  Comput.  Sci. 
3320  (2004)  184-187. 

[11]  J.R.  Ferguson,  J.M.  Fiard,  R.  Herbin,  Three-dimensional  numerical  simu¬ 
lation  for  various  geometries  of  solid  oxide  fuel  cells,  J.  Power  Sources  58 
(1996)  109-122. 

[12]  A.M.  Al-Qattan,  D.J.  Chmielewski,  S.  Al-Hallaj,  J.R.  Selman,  A  novel 
design  for  solid  oxide  fuel  cell  stacks,  Chem.  Eng.  Sci.  59  (2004)  131-137. 

[13]  R.  Saisset,  C.  Turpin,  S.  Astier,  B.  Lafage,  Study  of  thermal  imbalances  in 
arrangements  of  solid  oxide  fuel  cells  by  means  of  bond  graph  modeling, 
in:  Proceedings  of  the  IEEE  33rd  Annual  Power  Electronics  Specialists 
Conference,  vol.  1,  2002,  pp.  327-332. 

[14]  P.  Costamagna,  L.  Magistri,  A.F.  Massardo,  Design  and  part-load  perfor¬ 
mance  of  a  hybrid  system  based  on  a  solid  oxide  fuel  cell  reactor  and  a 
microgas  turbine,  J.  Power  Sources  96  (2001)  352-368. 

[15]  M.G.  Zhang,  X.G.  Wang,  W.H.  Li,  Nonlinear  system  identification  using 
least  squares  support  vector  machines,  in:  Proceedings  of  the  International 


H.-B.  Huo  et  al.  /  Journal  of  Power  Sources  162  (2006)  1220-1225 


1225 


Conference  on  Neural  Networks  and  Brain,  ICNN&B  05,  vol.  1  (13-15), 
2005,  pp.  414-418. 

[16]  Y.H.  Li,  S.S.  Choi,  S.  Rajakaruna,  An  analysis  of  the  control  and  operation 
of  a  solid  oxide  fuel-cell  power  plant  in  an  isolated  system,  IEEE  Trans. 
Energy  Convers.  20  (2)  (2005)  381-387. 

[17]  M.  Espinoza,  J.A.K.  Suykens,  B.  De  Moor,  Least  squares  support  vec¬ 
tor  machines  and  primal  space  estimation,  in:  Proceedings  of  the  42nd 


IEEE  Conference  on  Decision  and  Control  Maui,  Hawaii,  USA,  2003,  pp. 
3451-3456. 

[18]  C.-C.  Chang,  C.-J.  Lin,  LIBSVM:  a  library  for  support  vector 
machines,  2001.  Software  available  at  http://www.csie.ntu.edu.tw/~cjlin/ 
libsvm. 

[19]  T.  Hansen,  C.J.  Wang,  Support  vector-based  battery  state  of  charge  estima¬ 
tor,  J.  Power  Sources  141  (2005)  351-358. 


